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ABSTRACT 

We investigate the impact of feedback - from supernovae (SNe), active galactic nuclei (AGN) 
and a photo-ionizing background at high redshifts - on the neutral atomic hydrogen (HI) mass 
function, the 6j band luminosity function, and the spatial clustering of these galaxies at z=Q. 
We use a version of the semi-analytical galaxy formation model GALFORM that calculates 
self-consistently the amount of HI in a galaxy as a function of cosmic time and links its star 
formation rate to its mass of molecular hydrogen (H 2 ). We find that a systematic increase or 
decrease in the strength of SNe feedback leads to a systematic decrease or increase in the 
amplitudes of the luminosity and HI mass functions, but has little influence on their overall 
shapes. Varying the strength of AGN feedback influences only the numbers of the brightest 
or most HI massive galaxies, while the impact of varying the strength of photo-ionization 
feedback is restricted to changing the numbers of the faintest or least HI massive galaxies. 
Our results suggest that the HI mass function is a more sensitive probe of the consequences 
of cosmological reionization for galaxy formation than the luminosity function. We find that 
increasing the strength of any of the modes of feedback acts to weaken the clustering strength 
of galaxies, regardless of their Hi-richness. In contrast, weaker AGN feedback has little effect 
on the clustering strength whereas weaker SNe feedback increases the clustering strength of 
Hi-poor galaxies more strongly than Hi-rich galaxies. These results indicate that forthcoming 
HI surveys on next generation radio telescopes such as the Square Kilometre Array and its 
pathfinders will be exploited most fruitfully as part of multiwavelength survey campaigns. 

Key words: galaxies: formation - evolution - large-scale structure of the Universe - radio 
lines: galaxies 



1 INTRODUCTION 

Neutral hydrogen, both atomic (HI) and molecular (H2), plays a 
fundamental role in galaxy formation, principally as the raw mate- 
rial from which stars are made. The amount of neutral hydrogen in 
a galaxy at any given time reflects the complex interplay between 
processes that either deplete it, such as supernovae, or replenish it, 
such as gas cooling from hot atmospheres surrounding galaxies and 
by mergers with other gas-rich galaxies. By quantifying the prop- 
erties of HI in galaxies and noting how these properties vary with 
galaxy morphology and environment, we can glean insights into 
the physics of galaxy formation and test predictions of theoretical 
models. 

Most of what we know about the HI properties of galax- 
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ies comes from surveys of the nearby Universe ( z < 0.05) 
such as HIPASS (HI Parkes All-Sky Survey; see iMever et all 
20041) and ALFALFA (Arecibo Legacy Fast ALFA Survey; see 
Giovanelliet~al]|2005l) . Such sur veys have revealed t hat Hl-rich 
galax ies tend to be late-type (e.g. iKilborn et al .11 20021 ; lEvoli et all 
1201 lh, but it is comm on for early-type field galaxies to host HI 
(e.g. lserraet~al]|2012h ; that the HI mass function is well described 
by a Schechter function (e.g. lZwaan et alj2003l:lMartin et alj201C ) 
but its shape depends on environment (e.g. IZwaan et alj 12005 ; 
ISpringob et al]|2005l : IKilborn et alj|2009l) ; and that Hl-rich galax- 
ies are among the most weakly clustered galaxies known (e.g . 
IMever et alj |2007| ; iBasilakos etal] 120071 : IPassmoor etall 1201 lh . 
However, knowledge of the HI properties of galaxies will im- 
prove many-fold over the coming decade with the advent of 
next generation HI gal axy surveys on ASKAP (cf. Ijohnston et alj 
120081) , MeerKAT (cf. Ijonaj 120071) ) and ultimately the SKA it- 
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self (e.g. iBaugh et alj|2004 iBlake et alj|2004 IPower et al.ll20ld : 
lAbdalla et al.ll2010l : lKim et alfeoill) . ~ 

Future galaxy surveys on the SKA and its pathfinders are ex- 
pected to revolutionize our view of the HI Universe and so it is 
timely to ask precisely what these surveys can teach us about the 
physical processes that drive galaxy formation. In this paper we fo- 
cus on feedback - from stars in the form of winds driven by super- 
novae (SNe), from accreting super-massive black holes in the form 
of active galactic nuclei (AGN) heating, and from a photo-ionizing 
background in the early Universe - and evaluate how the global 
properties of HI in galaxies such as the HI mass function, which 
quantifies the number density of galaxies of a given HI mass, and 
the 2-point correlation funct ion and halo occupation distributions 
(HODs; e.g. lKimetal.l201 1), which quantify spatial clustering, are 
shaped by different sources of feedback. 

In galaxy formation models, the strengths of these processes 
and the interplay between them have been traditionally set by ex- 
amining the predictions of the galaxy luminosity function in the 
fej band, so we will also study how this statistic changes. To do 
this, we use the ve rsion of the semi -analytical galaxy formation 
model GALFORM o f ICole et all tod) as it has been developed by 
lLagos et"al] d2011al) : this calculates self-consistently the HI prop- 
erties of galaxies by splitting their interstellar media (ISM) into HI 
and H 2 pha ses using the empir ical relations of lBlitz & Rosolowskvl 
d2006h and lLerovetail J2008h . and links star formation in galax- 
ies not to their cold gas masses, as assumed in previous models 
(cf. lCole et alJl200d : lBa"ughll2006l), but to their H 2 m asses, as sug- 
gested by recent obs ervations (e.g. iBigiel et alj2"008h . As shown in 
lLagos et al] fcOllbl) . this model reproduces the observed HI mass 
function at z=0, accurately reproducing its amplitude and shape at 
intermediate and low HI masses. 

This is particularly interesting because we expect feedback to 
influence the global properties of galaxies selected either by their 
stellar mass or light. Previous studies have shown that SNe are 
pivotal in fixing the ampli t ude and slope of th e luminosity func- 
tion (e.g. ICole et al] |2000| : IBenson et alj|2003l) while AGN heat- 
ing suppresses the formati on of massive galaxies and dictates th e 
form of the bright end (cf. iBower eUdl2006l : ICroton et al.ll2006l) . 
In contrast, the influence o f SNe and AGN on the HIMF appears 
more subtle. For example, IPower et"al] fcOld) found that galaxy 
formation models that inc lude d or excluded AGN heating (such as 
|Pe Lucia & Blaizoll2007l and lBaugh et alj|2005l respectively) pre- 
dict HI mass functions that reproduce equally well the observed 
z ~ data. The HI mass function offers the possibility of placing 
stronger constraints on the strength of feedback in semi-analytical 
models, which previously u s ed only the optical properties of galax- 
ies (e.g. IBenson et al.ll2003l ; lB~ower et alj|201(]|) . At the same time, 
the variation of clustering strength with galaxy properties provides 
us with important clues about the physics of galaxy formation. Any 
discrepancy between observational measurements of clustering and 
the predictions of galaxy formation models indicates the need to 
improve the models, either by refining the modelling of the phys- 
ical processes included or by considerin g the addition of further 
effects (see, for example. iKim et all2009t) . 

The structure o f the paper is as foll ows. In Sj2] we provide a 
brief overview of the lLagos et al .1 d201 lat) model and describe how 
the different modes of feedback (SNe, AGN and photo-ionization) 
are implemented in GALFORM. In ij3]we show how the different 
forms of feedback influence both the galaxy luminosity function in 
bj band and the HI mass function, and we deduce using likelihood 
maximization precisely what information we glean from them. In 
Sj4] we investigate how the spatial clustering of Hi-poor and HI- 



rich galaxies are influenced by feedback by inspecting the 2-point 
correlation function and HODs. Finally, in Sj5] we summarize our 
results and discuss their implications for testing the modelling of 
feedback with forthcoming HI surveys. 



2 MODELLING 

In this section we first briefly give an overview of the version of 
GALFORM used in this paper (§ 12. It , Because the focus of our pa- 
per is on the impact of feedback on the large-scale HI properties 
of galaxies, we summarize the modes of feedback that are imple- 
mented in GALFORM - from SNe (§[2T2j, AGN (§|23) and photo- 
ionization after cosmological reionization (§ |2.4| l. 



2.1 The GALFORM galaxy formation model 

We use t he version of GALFO RM (cf. ICole et alj|2000h that is de- 
scribed in lLagos et alj J201 1 al Tbh to predict the properties of galax- 
ies forming and evolving in t he ACDM cosmology adopted for 
the Millennium Simulation (cf. ISpringel et al. ll2005li T l GALFORM 
models the key physical processes of galaxy formation, includ- 
ing the gravitationally driven assembly of dark matter halos, radia- 
tive cooling of gas and its collapse to form centrifugally supported 
discs, star formation, and feedback from supernovae (SNe) and ac- 
tive galactic nucl e i (AGN ). 

lLagos et al] d2011al) extended GALFORM by modelling the 
splitting of cold gas in the interstellar medium (ISM) into its HI 
and H2 components and by linking e xplicitly star format ion to 
the amount of H2 present in a galaxy. lLagos et al.1 J201 ltah com- 
pared the model predictions obtained using empirically and theo- 
retically derived star f ormation laws (cf. lBlitz & Rosolowsky|20o3 : 
iKrumholz etai]|2009l) with a variety of observations (e.g. the HI 
mass function, 12 CO (1-0) luminosity function, and correlations 
between the ratio H2/HI and s tellar and cold gas masses ) and 
fou nd that the empiri cal law of iBlitz & Rosolowskvl J2006h (see 
also lLerov et al.|[2008h is favoured by these data. This law is of the 
form 

SsFR = ^SF fmol Sgaa, (1) 

where Esfr and E gas are the surface densities of the star formation 
rate (SFR) and total cold gas mass respectively, ^sf is the inverse 
of the SF timescale for the molecular gas and f mo i = E mo i/E ga s 
is the molecular to total gas mass surf ace density ratio. Imp ortantly 
for the work we present in this pap er, lLagos et all feOllbl) showed 
that the lBlitz & Rosolowskvl | |2006|) law is able to reproduce the HI 
mass function at z=Q at intermediate and low HI masses. This is be- 
cause it suppresses star formation in lower mass galaxies, thereby 
reducing SNe feedback and allowin g these galaxies to re tain larger 
gas reservoirs. Note that we use the lLagos et al.1 ( l20Tlbl) as the de- 
fault model in this paper. 



2.2 Supernovae (SNe) Feedback 

SNe act to reheat cold gas and eject it from the galaxy disc at a rate 



1 Recall that the cosmological parameters adopted for the Millennium 
Simulation are the total matter density = 0.25, the baryon density 
f!b = 0.045, the vacuum energy density Q\ = 0.75, the Hubble pa- 
rameter Hq = 100ft kms - 1 Mpc -1 with /i=0.73, the primordial scalar 
spectral index n B = 1 and the fluctuation amplitude erg = 0.9. 
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(2) 



where ip is the instantaneous star formation rate and /3 regulates 
feedback efficiency, 



P = (V disk /t/ hot )- Qh 



(3) 



Here Vhot and «hot are adjustab le parameters th at govern the 
strength of the SNe feedback (cf. ICole et alj|200(jl) and V dis k is 
the circular velocity of the galactic disc; the values in the default 
ILagos et"al] d2011bl) model are Vhot=485 km/s and ahot=3.2. We 
expect SNe feedback to be most effective in low-mass dark matter 
halos because gas is expelled efficiently from their shallow poten- 
tial wells. From Eq[3] the mass loading of the SNe wind is largest 
in such halos and this suppresses the formation of faint, low-mass 
galaxies and ensures that the faint end of the galaxy luminosity 
function is recovered b y GALFORM, in good agreement with obs er- 
vational estimates (e.g. lBlanton et al ]|200ll ; lNorberg et al.l2002h . 



2.3 AGN feedback 

SNe feedback is less effective in large circular velocity halos (ac- 
cording to Eq[3]l. The gas ejected from low mass halos can poten- 
tially cool and turn into stars in more massive halos, generally lead- 
ing to the formation of too many massive galaxies. In the absence of 
strong feedback in massive dark matter halos, GALFORM predicts 
an overabundance of bright galaxies and a gal axy luminosity func - 
tion whose bri ght end slope is too steep (cf. iBenson et al. I l2003l) . 
For this reason IBower et al.l d2006l) extended GALFORM to include 
AGN feedback, which suppresses cooling flows in massive halos. 
Physically this occurs because AGN are fuelled by super-massive 
black holes that grow by luminous accretion. 

In GALFORM AGN feedback is modelled by assuming that 
the hot gaseous atmospheres surrounding massive galaxies are in 
quasi-hydrostatic equilibrium when the cooling time at the cooling 
radius, i C ooi(Vcooi), exceeds a multiple of the free-fall time at the 
cooling radius, tfi(r coo i), i.e. 

1 



^cool (^"cool ) ^* 



•^cool 



-£ff (Vcooi); 



(4) 



where a CO oi is a free parameter that regulates the strength of AGN 
feedback.! increasing a coo i means that AGN feedback becomes 
effective in lower mass dark matter halos. We adopt a v alue of 
a cool =0.58, which is the default in the lLagos et al.l O01 lbh model. 
We define the cooling time at halocentric radius r as 



f ^ 3 



firripkBTg, 



2p gas (r)A(T gas ,Z gas ) 



(5) 



where p gas , T sas and Z gaa are the gas density, temperature and 
metallicity respectively, A is the cooling functiorQ ft is the mean 
molecular weight, m p the proton mass and ke is Boltzmann's con- 
stant. We define the free-fall time at halocentric radius r by 



2GM(x) 



Ax 



-1/2 



dy. 



(6) 



where M(r) is the total mass (baryonic and dark matter) interior 
to r. The cooling radius, r CO oi> is the radius at which t coo i matches 
the age of the halo, which is obtained from the halo merger history. 

We note that the model predicts that generally black holes gain 
most of their mass by accretion of gas driven to the centre of the 



2 We note that Eq 2 of lBower et alj feOOd) contains a typographical erro r. 

3 We use the tabulated cooling functions of Sutherlan d & Dopit3 J 1993b - 



galaxy via disc instabilities (cf. iFanidakis et all 1201 ll) . Accretion 
of gas that cools from the hot atmosphere is relatively unimportant 
until late times, when the growth rate of the black hole scales as 
icooi / c 2 , where L coo i is the luminosity of the cooling gas and c is 
the speed of light. 



2.4 Photo-ionization feedback 

Photo-ionization is predicted to have a dramatic impact on star for- 
mation in low-mass galaxies. This is because the presence of a 
photo-ionizing background both modifies the net cooling rate of 
gas in halos by removing the "hy drogen peak" in the cooling curve 
(cf. Fig. 1 of IBenson et al. I l2002l) and increases the temperature of 
the intergalactic medium (IGM) such that its thermal pressure pre- 
vents gravitational collapse of baryon s onto low-mass halos (e.g. 
lEfstathioull992l ; IOkamoto et alj2008l) . As a result, only those low- 
mass halo s that contained c old gas prior to re-ioniza tion can form 
stars ( e.g. lrloeft et alfeoOot) . GALFORM includes the IBenson et alj 
j2002l) prescription for suppressing the cooling of halo gas onto the 
galaxy - this occurs if the host halo's circular velocity Vdic lies be- 
low a threshold V cu t at redshift z < z cu t- The values in the default 
ILagos etal.ld201 lbh model are V cu t=30 km/s and z cut = 10. The 
default Vcut is in g ood agreement wit h th e results of hydrodynam - 
ical simulations bv lHoeft et ail d2006h and lOkamoto et al 1 12001) . 



3 THE GALAXY LUMINOSITY FUNCTION AND HI 
MASS FUNCTION 

In this section we investigate how the galaxy luminosity function 
in bj band and the HI mass function at z=0 are influenced by 
feedback from supernovae (SNe), active galactic nuclei (AGN) and 
a photo-ionizing background. We use Monte Ca rlo halo merger 
trees g enerated using the algorithm described in IParkinson et alj 
d2008l) to predict galaxy properties. The Monte Carlo trees are 
tuned to match A^-body merger trees extracted from the Millen- 
nium Simulation. Note that we express HI masses in units of 
h~ 2 M@ rather than ft _1 M0, the mass unit used in simulations. 
This ensures that the observational units (which depend upon the 
square of the luminosity distance) are matched, but it introduces 
an explicit dependence on the dimensionless Hubble parameter, 
for which we assume the value of fe=0.73 used in the Millenn ium 
simulation. We refer the reader to § 2. 1 o f lLagos et alj fcOllbh for 
further discussion of this approach. 

Intuitively, we expect SNe feedback to be particularly damag- 
ing in low-mass galaxies because gas is expelled easily from their 
shallow potential wells and so star formation is shut down. In con- 
trast, SNe feedback is ineffective in high-mass galaxies, which re- 
side at the centres of massive dark matter halos with deep poten- 
tial wells; here a far more energetic source of feedback, namely 
AGN heating, is required, which acts to suppress the formation 
of very massive bright galaxies. The photo-ionizing background 
suppresses the collapse of gas onto low mass dark matter halos, 
thereby preventing galaxies from forming. The influence of these 
distinct processes on the galaxy luminosity function is straightfor- 
ward to predict - photo-ionization and SNe feedback will govern 
the amplitude and shape of the faint end, while AGN heating dic- 
tates the amplitude and shape of the bright end. Predicting their 
i mpact on the HI mass function is less straightforward - as noted 
in ILagos etaljfcfjl lbh . the most Hl-rich galaxies are likely to be 
low-mass systems, in which case the imprint of SNe feedback may 
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span the mass function, whereas the influence of AGN heating may 
be negligible. 

We investigate these ideas in FigQ] in which we explore the 
impact of SNe driven winds (upper panels), AGN heating (mid- 
dle panels) and a photo-ionizing background (lower panels) on the 
galaxy luminosity function in bj band (left panels) and HI mass 
function (right panels) at z=Q. For comparison we plot the bj band 
galaxy luminosity function (open triangles) estimat ed from the 2 
degree Field Galaxy Redshift Survey (2dFGRS) by iNorberg et al.l 
d2002l) in the left panels, and the HI mass functions deduced from 
the HI PASS (open triangle s) a nd ALFALFA (fille d squares) sur- 
veys bv lZwaanetai] d2005h and lMartin etal] fcOlCh respectively. 

Impact of SNe feedback: We vary the parameter Vh ot , which gov- 
erns the mass ejection rate, M e ject, assuming a fixed value of a no t, 
between 300km/s (less effective) and 700km/s (more effective); 
recall its default value is 485 km/s. By varying Vhot we are in 
essence changing the mass of the system in which SNe feedback 
is effective. On the other hand varying otbot would result in the 
same galaxies being affected by SNe, but with different mass load- 
ing in the winds of reheated material. The upper panels indicate that 
varying Vhot in this way leads to systematic and measurable offsets 
in both the luminosity function and HI mass function - stronger 
(weaker) SNe feedback decreases (increases) the number of galax- 
ies of a given luminosity (by ~ 0.5-1 dex) or HI mass (by ~ 0.5 
dex). Interestingly, weaker feedback leaves the shape of the lumi- 
nosity and HI mass functions broadly unaffected, whereas stronger 
feedback leads to slightly steeper slopes at the faint/low-mass and 
bright/high mass ends. Red dashed and blue dotted curves cor- 
respond to models with strong (600km/s) and weak (300km/s) 
feedback; we use these as reference models in §f4] where we exam- 
ine the clustering of Hi-rich and Hi-poor galaxies. 

Impact of AGN feedback: Here we vary the parameter q coo i be- 
tween 0.4 ^ a coo i ^ 0.8; recall its fiducial value is a CO oi = 0.58. 
Increasing a CO oi decreases the halo mass at which AGN heating 
starts to become effective. The middle panels of FigQ]demonstrate 
that AGN heating has little influence beyond the bright end of the 
galaxy luminosity function or the high-mass end of the HI mass 
function, regulating both the sharpness of the transition between 
faint/bright and low/high mass galaxies and the position at which 
it occurs. The number of bright galaxies is very sensitive to the 
strength of AGN feedback, varying by as much as a factor of ~ 100 
for the most luminous systems. Interestingly AGN feedback also 
affects the numbers of galaxies with high HI masses, but the effect 
is only apparent in the cases of weakest feedback. Red dashed and 
blue dotted curves correspond to models with strong (a coo i = 0.8) 
and weak (a coo i = 0.4) AGN feedback, which, as before, we use 
as reference models in §f4] 

Impact of Photo-ionization: Two parameters, Vut and z cu t, gov- 
ern photo-ionization in the model. We find that z cu t, the redshift at 
which reionization occurs, has little influence on either the galaxy 
luminosity function or HI mass function at z=0 - varying z cut be- 
tween 6 ^ z cut ^ 12, the range suggested by measurements of 
polarisation in the cosmic microwave backgro und radiation and the 
Gunn-Peterson effect in distant quasars (e.g. iGreiner et al7ll2009l ; 
iKomatsu et al. ll2oTll) . produce luminosity and HI mass functions 
that are indistinguishable, to within the line width. For this rea- 
son we concentrate on V cu t, which we vary between 30km/s ^ 
Vcut ^ 90km/s. The minimum progenitor mass we consider is 
5 x 10 s /i -1 © to explore low V cu t values (see also HI mass func- 



tion in figure. [T]which has lower HI mass limit than other cases of 
HI mass function, indicated by the vertical dotted line). Increas- 
ing Vcut suppresses cooling in all halos with circular velocities 
Vdr ^ Vcut at z < Zcut and its impact is most pronounced at 
the faint end of the galaxy luminosity function and the low-mass 
end of the HI mass function. The red dashed curve corresponds to a 
model with strong (50km/s) photoionization feedback, which we 
use as a reference model in §|4] 

Constraints on Galaxy Formation Models: What information 
can we glean from the galaxy luminosity and HI mass functions 
about the various modes of fee dback? We investigat e this using 
the galaxy lumin osity function of.Norberg et al. ( 2002 ) and the HI 
mass function of IZwaan et alj l l2005t) . Although the [Martin et al .1 
ALFALFA mass function ext ends to lower HI mas ses, it has 
more limited sky coverage than the IZwaan et alj J2005h HIPASS 
mass function and so is more susceptible to the effects of cosmic 
variance. We quantify the constraints offered by the observational 
data by calculating the ratio of likelihoods L/Ldefauit for each 
models with respect to default model, 

L = S|log(odi)-log(mdi)| 

Ldcfauit S| log(odi) - log(mdi)| default ' 

Here we sum over the difference in log values of the observational 
data (odi) and the predictions (mdi) for each model and compare 
with the sum over this difference for the default model. Our re- 
sults are plotted in Fig [2] where solid (dotted) curves indicate how 
L/Ldeiau\t varies with SNe, AGN and photo-ionization feedback 
strength (left, middle and right panels) for the galaxy luminosity 
function (HI mass function). 

Formally, the models give a poor fit to the observational data. 
There are several reasons for this. Firstly, the parameters we vary 
are physical, not statistical, and affect the predictions of the model 
in a complex way, due to the interplay between different physical 
processes. With a parametric functional form, such as the Schechter 
function, the consequences of changing a parameter are transpar- 
ent. Secondly, GALFORM predicts a whole range of galaxy prop- 
erties, and the luminosity function or HI mass function represent 
just one way of describing the population of model galaxies. The 
parameters of the default model were chosen as a compromise in 
reproducing a variety of observational data, not just the luminosity 
function or mass function. If we chose to fit the data considered 
here in isolation, with complete disregard for any other dataset, it 
would be possible to improve the fits. Nevertheless, it is not clear 
that we would obtain a formally acceptable fit. Again, the reason 
for this is that GALFORM is a physical model rather than a para- 
metric fitting formula. The shortcomings of the model could be en- 
cod ed in a "model disc repancy" that broadens the errors in a \ 2 fit 
(see lBower et alj2010l for a more formal approach to this). Finally, 
the observational errors do not include all of the contributions (e.g. 
sample variance) and exclude systematics (e.g. the k-correction in 
the case of the luminosity function). The key point we aim to get 
across from this plot is the shape of the L/Ldcfauit curves: How 
rapidly do we move away from the quality of the fit in the default 
model on varying the parameter? If the value of L/Ldefault varies 
rapidly on comparing the model predictions to a particular dataset, 
then that dataset could in principle provide a tighter constraint on 
the parameter value. 

We find that the galaxy luminosity function shows a strong 
dependence on the strength of SNe feedback and a moderate de- 
pendence on the strength of AGN feedback, with a narrow range 
of preferred values close to their default values; the HI mass func- 
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Figure 1. Impact of feedback from SNe, AGN and a photo-ionizing background (top, middle and bottom rows) on the predicted b j -band galaxy luminosity 
function (left pa nels) and HI mass fu nction (right panels); symbols correspond to data from the 2dF Galaxy Redshift Sur vey (open triangles ; cf. Figure 1 1 
and Table 1 of iNorberg et afll2002l> . HIPASS (open triangles; cf. IZwaan et alj|2005h and ALFALFA (filled squares; cf. iMartin et alj|201(j) . SNe strength 
varies between 300km/s < Vhot ^ 700km/s, assuming fixed ai lot ; AGN strength between 0.4 ^ a CO ol ^ 0-8; and photo-ionization strength between 
30km/s ^ V C ut ^ 90km/s. The darker the line, the closer the strength of the feedback parameter to its default value. Red dashed and blue dotted curves 
correspond to models with strong and weak feedback; the dotted vertical line in the lower right panel indicates the lower HI mass limit in the galaxy formation 
model we use them as reference models in §|4] 
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tion offers a weak constraint on SNe feedback and no constraint 
on AGN feedback. In contrast the HI mass function is sensitive to 
the strength of the photo-ionizing background, and indicates that 
values of V cu t larger than ~ 70 km/s can be ruled out. 

Smaller values of Kut appear to have the same likelihood, 
but this can be understood as an artifact of the effective resolu- 
tion of our merger trees. We use Monte Carlo merger trees to ex- 
plore the impact of varying Kut, as this allows us to vary the mass 
resolution of the progenitor halos to ensure that we can resolve 
the halos which would be affected by the choice of this parame- 
ter. The minimum progenitor mass we consider is 5 x 
We can obtain tighter constraints on the minimum value of V cut 
by using Monte Carlo trees with higher effective resolution, but 
we defer this to a forthcoming paper. Our results suggest that the 
HI mass function, measured at 2=0, can provide an important con- 
straint on the minimum halo mass that contributed to reionization, 
and complements estimates estimates of t he minimum halo mass 
based on high redshift observations (e.g. Choudhurv et al .1 1200 8l ; 
ISrbinovskv & WvithtfeOloHMufioz & Loebll201ll) 



4 SPATIAL CLUSTERING OF HI IN GALAXIES 

We now examine how feedback influences the spatial distribution 
of HI galaxies, which we quantify using the 2-point correlation 
function and the ha l o occupation distributio n (hereafter HOD; cf. 
Benson et alj |200C| : IPeacock & Smith! l2000l : iBerlind & Weinberd 



20021) . The clustering predictions for fej-band selected sample 



were studied in Kim et al. (2009). The HOD gives the mean num- 
ber of galaxies as a function of halo mass, assuming pre-defined 
selection criteria (e.g. colour, cold gas mass), and can be sepa- 
rated into contributions from central galaxies and their satellites. 
It is an output from the semi-analytical model. HODs for opti- 
cally selected samples of galaxies have been studied extensively 
and are typically parameterized as step functions, which represent 
central ga laxies, and power-law components, which represent s atel- 
lites (e.g. IPeacock & Smith! [200ol ; IBerlind & Weinberg! |2002|) . In 
contrast, relatively few studies have been made of HODs for HI se- 
lected samples of galax ies (e.g. lWvithe & Brownl20ld : lMarrn et al.1 
l2010l :l Kim et al]|201ll) but they indicate that parameterizations of 
the kind adopted for optically select ed galaxy sample s need to be 
modified for HI selected samples (cf. lKim et al J201 lh . 

Note that we use iV-body halo merger trees drawn from the 
Millennium Simulation to obtain the spatial information required 
to study clustering. The resolution limit of the Millennium Simula- 
tion imposes a limit on the halo mass that we can reliably resolve 
of 10 10 h' 1 M Q , which in turn sets a limit of Mm ~ 10 8 5 /i _2 M o 
on the HI mass that we can track. 

In Fig [3] we explore the impact of SNe driven winds (up- 
per panels), AGN heating (middle panels) and photo-ionization 
(lower panels) on the 2-point correlation function (left panels) and 
the HOD (right panels) at z=Q. In the top part of each panel we 
plot the relevant statistic for galaxies selected according to their 
high or low HI mass, while in the bottom part we plot the ratio 
of the statistic with respect to its value in the default model for 
galaxies of the same HI richness. Here we define Hi-rich galax- 
ies to have Aim > 10 9 ' 25 h~ 2 M@ while Hi-poor galaxies have 
lO 8 ' 5 /i," 2 M < Mm < 10 9 ' 25 /i" 2 M Q . For clarity, we concen- 
trate on the typical weak and strong feedback models (blue and 
red curves respectively) that we considered in Fig[T] This is neces- 
sary because we expect environmental effects to become apparent 



when we measure clustering properties, and we wish to disentan- 
gle the effects of environment from feedback. We show also the 
2-point correlation function of gas ri ch galaxies at z=0 m easured 
from the HIPASS survey, presented in iMever et all d2007h . and the 
predictions for the 2-point correlation function and the HOD in the 
default model, which is in good agreement with the HIPASS result. 

Impact of SNe feedback: The default model predicts that Hi-rich 
galaxies (black solid curve) should be more strongly clustered than 
Hi-poor galaxies (red dotted curve), which we show in the top left 
panel. Here we see also that strong SNe feedback has negligible 
impact on the clustering strength at large separations of both HI- 
rich and Hi-poor galaxies (dotted-dashed and dashed curves re- 
spectively), but it influences the clustering strength of both Hi-rich 
and Hi-poor galaxies at separations < lh^ 1 Mpc (short dashed 
curve). Reducing the strength of SNe feedback leads both Hi-poor 
and Hi-rich galaxies to cluster more strongly at small separations 
- < l(3)h~ Mpc for Hi-rich (Hi-poor) galaxies. Fig. [3] shows 
that the clustering amplitude of both Hi-poor and Hi-rich samples 
becomes similar - as can be most easily seen in the lower panel, 
which shows the ratio of the clustering with respect to that of the 
corresponding HI sample in the defa ult model. This tr end can be 
understood by inspecting Figure 3 of iKim et alj feoill) . in which 
the cold gas mass in galaxies is plotted against their host dark mat- 
ter halo mass. Beyond the mass at which AGN heating stops cool- 
ing flows, there is a sudden break in the cold gas mass - halo mass 
relation and a dramatic increase in the scatter. Above this break in 
mass, the presence of gas is due to accretion of halos of lower mass, 
which have their own gas reservoirs, onto the main halo. By reduc- 
ing the strength of SNe feedback, more gas is able to stay in the 
cold phase in these progenitors and gas is depleted more slowly in 
the current host. Because there is a large scatter in galaxy cold gas 
mass for a given halo mass, there is a shift in the typical host mass 
for both Hi-poor and Hi-rich galaxies, leading to an increase in the 
clustering amplitude. 

In the top right upper panel we show the predicted HODs as 
a function of HI richness (dashed and dotted curves) for all galax- 
ies. Thick (thin) curves indicate strong (weak) SNe feedback cases 
for Hi-poor and Hi-rich galaxies. As we would expect, stronger 
SNe feedback leads to fewer galaxies when compared to the weaker 
feedback case, spanning the full range of halo masses and regard- 
less of the HI richness of the galaxies. In the top right lower panel 
we plot the ratio of the number of galaxies for the strong and weak 
SNe feedback cases compared to the default model. Again this re- 
veals that SNe feedback affects the numbers of galaxies regardless 
of host halo mass, but interestingly it shows that the influence of 
SNe feedback is greater on the numbers of Hi-rich galaxies com- 
pared to Hi-poor galaxies. 

Impact of AGN feedback: In the middle left panel we show that 
the clustering of galaxies selected by their HI richness is largely 
insensitive to the strength of AGN feedback. If AGN feedback is 
weak, its effect on the clustering strength of Hi-poor galaxies is 
negligible and enhances the clustering of Hi-rich galaxies only at 
small separations. If the AGN feedback is strong, it suppresses the 
clustering of both Hi-rich and Hi-poor galaxies to a similar extent. 

As before, we show in the middle right upper and lower panels 
the predicted HODs as a function of HI richness (dashed and dot- 
ted curves) for all galaxies, where strong (weak) AGN feedback is 
indicated by thick (thin) curves. The lower panel shows that AGN 
feedback affects halos with masses in excess of lO'^r'Mo and 
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Figure 2. Values of i/idcfault obtained for the SNe, AGN and photo-ionization feedback models (left, middle and right panels), where i/idefault is 
calculated using equation [7] Solid (dotted) curves show how L / Ldefault for the galaxy luminosity function in b j band (HI mass function) as we vary 
systematically the SNe, AGN and photo-ionization strength; filled squares (circles) indicate the values for the models that have been run. For reference we 
show also a dashed horizontal lines corresponding to £/£default=l- See text for further details. 



that it has a greater impact on the numbers of Hi-rich galaxies than 
Hi-poor galaxies. 

Impact of Photo-ionization: Increasing the strength of photo- 
ionization has negligible impact on the clustering strength of HI- 
rich galaxies, but leads to an increase in the clustering strength of 
Hi-poor galaxies. This is because increasing the strength of photo- 
ionization corresponds to an increase in the circular velocity scale 
Vcut of dark matter halos in which galaxies can form. 

We show predicted HODs as a function of HI richness (dashed 
and dotted curves) for all galaxies in the bottom right upper and 
lower panels. Here we consider only the case of strong photo- 
ionization feedback because the typical halo mass affected by weak 
photo-ionization feedback (V^ ut =30 km/s) drops below the resolu- 
tion limit of the Millennium Simulation at early times. Thick lines 
indicate strong photo-ionization feedback whereas thin lines show 
the default model for Hi-poor and Hi-rich galaxies. It is readily 
apparent that photo-ionization feedback has little influence over 
the HODs of Hi-rich galaxies - the strong feedback case is in- 
distinguishable from the default model. The lower panel shows 
that photo-ionization feedback affects only the lowest mass halos 
(lO lo - 7 /i _1 M )at2=O. 



5 SUMMARY 

We have used the ver sion of the GALF ORM semi-analytical galaxy 
formation model (c f. ICole et al] |200Cf) that has been extended by 
lLagos et"al] J201 lal) to study how feedback from SNe, AGN and 
photo-ionization shapes both the galaxy luminosity and HI mass 
functi ons, and the spatial clustering of HI galaxies. The advantage 
of the lLagos et alj fcOllal) model is that it is one of the first to pre- 
dict explicitly and in a self-consistent ma nner how much HI should 
reside in a galaxy at a given time (see also lFu et al.l20 10T) and it has 
reproduced successfully the observed HI ma ss function at low-to- 
intermediate masses (cf. lLagos etai]|2011bh . The main results of 
our study can be summarised as follows; 

• Feedback from SNe regulates the amplitude and shape of both the 
galaxy luminosity function and the HI mass function. The effect 
is systematic, with more (less) efficient SNe feedback leading to 



fewer (more) galaxies of a given luminosity or HI mass. It also 
regulates the amplitude of the 2-point correlation function of HI- 
rich galaxies at small separations and the numbers of galaxies per 
halo mass as quantified by the HOD, regardless of their HI richness. 

• Feedback from AGN has little effect on either the faint-to- 
intermediate luminosity end of the galaxy luminosity function or 
low-to-intermediate mass end of the HI mass function - its impact 
is stronge st for the brightest galaxies, as previo us studies have ar- 
gued (e.g. lCroton et al]2006l ; lBower et ai1l2006l) . but it also affects 
the number of highest HI mass g alaxies. A similar result has been 
reported by iFabello et alj J201 it) , and it also helps to explain the 
cold g as mass -halo mass relation presented in Fig. 3 of iKim et alj 
d20 1 lh for the lBower etal] feOOfj) model. In contrast, its impact on 
the clustering of HI galaxies is minor or negligible. 

• Feedback from photo-ionization is most pronounced for the 
faintest and most Hi-poor galaxies. The redshift at which reion- 
ization occurs is not important - but the mass scale, as measured 
by Kut, is, and its influence is evident in both the shape and am- 
plitude of the faint end of luminosity function and the low-mass 
end of the HI mass function. Our analysis suggests that the low-to- 
intermediate mass end of HI mass function offers the potential to 
constrain the models of photo-ionization. Interestingly, we find that 
the locally measured HI mass function can constrain the minimum 
mass of dark matter halos that could have hosted galaxies that con- 
tributed to reionization. This complements existing studies of the 
high redshift Universe that have used, for exampl e, the electron 
scattering optic al depth (e.g. Choudhurv et al.|[2008l) . the Lyman-a 
forest and (e.g lSrbinovskv & Wvithe 20101) and the UV luminos- 
ity function (e.g. Munoz & Loebl 20 1 1 ) to estimate the minimum 
halo mass. Our analysis suggests that the circular velocities of ha- 
los in which galaxy formation was suppressed cannot be larger than 
~70 km/s. 

• Strong modes of feedback, as they are implemented in GALFORM, 
act to suppress the clustering strength of galaxies, regardless of 
their Hi-richness. Hi-poor galaxies cluster more strongly than HI- 
rich galaxies if the strength of SNe feedback is weakened, whereas 
there is little change in the clustering of either Hi-rich or Hi-poor 
galaxies if the strength of AGN feedback is weakened. 

Our study suggests that forthcoming HI galaxy surveys can be most 
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Figure 3. The impact of SNe, AGN and photo-ionization feedback (top, middle and bottom rows) on the 2-point correlation function (left panels) and the 
halo occupation distribution (right panels). In the upper part of each panel we plot the relevant statistic for Hi-poor and Hi-rich galaxies, while in the lower 
part we plot the ratio of the Hi-poor or Hi-rich statistic with respect to its value in the default (Hi-poor or Hi-rich) model. Here Hi-rich galaxies have masses 
Mhi > 10 925 Mq and Hi-poor galaxie s have HI masses between 1O 8,5 M0 and 1O 925 M0. Filled triangles show the 2-point correlation function of gas-rich 
galaxies obtained bv lMever et alj J2007I) . and assorted line types refer to the different models, as indicated by the legend in the left panels. See text for further 
details. 
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fruitfully exploited scientifically if they are part of a multiwave- 
length campaign. Such campaigns are being planned; for example, 
DI NGO (Deep I nvestigation of Neutral Gas Origins) on ASKAP 
("cf . lMeverll2009l) will probe the galaxy population out to z < 0.4 
and will be combined the GA MA (Galaxy And Mass Assembly) 
survey (cf. |Priveretal.ll20T J). These data-sets will allow us to 
study not only the HI properties of galaxies that host AGN or young 
starbursts, but also the effects of environment and evolution with 
redshift. 
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